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Abstract: The final state of Au + Au collisions at \/s = 130 AGeV at RHIC has been 
reconstructed within the framework of the Buda-Lund hydro model, by performing a 
simultaneous fit to preliminary PHENIX and STAR data on two-particle Bose-Einstein 
correlations and identified single particle spectra, and the Hubble constant is determined 
to be H = {ut) = 0.77 ±0.09. 


1. Introduction 

The reconstruction of hadronic final state from the measured single particle spectra and 
two-particle correlation functions is of great current research interest in high energy heavy 
ion collisions, in order to identify one or more new phases of hot and dense hadronic 
matter in the collisions of the biggest available nuclei at the highest available bombarding 
energies. It has been expected [||, |^, that for a Bjorken type of initial condition and 
for a long-lived, soft, transient Quark Gluon Plasma phase, pions are evaporated from a 
predominantly longitudinally expanding, transversely almost expansionless fireball during 
a long period of time and this signature of the QGP phase can be observed experimentally 
from the analysis of two-pion Bose-Einstein correlation functions. 

An alternative approach is to reconstruct the hadronic final state from the measureable 
single-particle spectra and two-particle correlation functions. From this reconstructed final 
state and the knowledge of the equation of state of hot and dense hadronic matter (e.g. 
from lattice QCD calculations) one can, in principle, reconstruct the initial state of the 
reaction by running the (relativistic) hydrodynamical equations backwards in time, and 
determine if this initial state had been in the QGP phase or not. Here we report on such 
a reconstruction within the framewo rk of the Buda-Lund hydro model. 
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2. Buda-Lund hydrodynamic model 


The Buda-Lund hydro parameterization (BL-H) has been developed in refs. |Q, ^ to de¬ 
scribe particle correlations and spectra in central heavy ion collisions at CERN SPS. 

We assume that high energy heavy ion collisions can be interpreted within the core¬ 
halo picture Q where the center of the interactions is the core which is surrounded by 
the decay products of long lived resonances, corresponding to the halo. The Buda-Lund 
hydro parameterization describes the core of the particle emitting region as a cylindrically 
symmetric, finite hydrodynamically expanding system |Q, with emission function 


Sc{x,p) (Tx = 


(2vr)' 


exp 


u^{x)p^ Pl{x) 


( 2 . 1 ) 




T{x) T{x) 

where the subscript c refers to the core of collision (surrounded by a halo of long lived 
resonances). The degeneracy factor is denoted by g, the four-velocity field is denoted by 
u^{x), the temperature field is denoted by T{x), the chemical potential distribution by 
p,{x) and s = 0, — 1 or 1 for Boltzmann, Bose-Einstein or Fermi-Dirac statistics. 

The particle flux over the freeze-out layers is given by a generalized Cooper-Frye factor, 
assuming that the freeze-out hypersurface depends parametrically on the freeze-out time 
r and that the probability to freeze-out at a certain value is proportional to H{t), 

pPd‘^Yi^{x) = mtcosh.{ri — y]H{T)dT T^dr] dcx dry, ( 2 . 2 ) 

H{t) = (27rAr^)“^/^exp [-(r - ro)^/(2Ar2)] , (2.3) 

where we keep only the mean and the variance of the effective prop er-time dist ribution 
H{t). The transverse mass and coordinate are m* = ^m? + p%+ Py and rt = 
the rapidity y and the space-time rapidity rj are defined as y = 0.5 log [{E + pz)/{E — pz)] 
and rj = 0.5 log [{t + z)/{t — z)]. 

The distributions of 1/T(x) and pi{x)/T{x) are parameterized within BL-H as 
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(2.5) 


The central temperature and chemical potential at the mean freeze-out time are denoted 
by To = T{rx = = 0;t = tq) and = 0 ;t = tq). With the surface 

temperature T^ = T{rx = fy = Rg,x = tq) and the temperature after freeze-out, Tt = 
T{rx = Cy = = tq + \/2At), the relative transverse and temporal temperature decrease 

are introduced, see refs. 0, H, for further details. The variation of the chemical 

potential in coordinate space is related to the finiteness of the density profile in the core. 

The four-velocity u^{x) of the expanding matter is assumed to have the form|Q, |l^, ]R|; 


u^{x) = |^cosh[77] cosh[? 7 t], sinh[r/t] —, sinh[r/t] —, sinh[r?] cosh[r 7 t 
sinh[T/t] = {ut)rt/RG, 


( 2 . 6 ) 
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where Rq stands for the transverse geometrical radius of the source. Such a flow profile, 
with a time-dependent radius parameter Rq, was recently shown to be an exact solution of 


relativistic hydrodynamics of a perfect fluid at a vanishing speed of sound |14|. It turned 
out 1^, , that the flow field is a generalized Hubble flow and the average transverse flow 

at the geometrical radius is formally similar to Hubble’s constant that characterizes the 
rate of expansion in our Universe, {ut) = Rg = H. This emphasizes the similarity between 
the Big Bang of our Universe and the Little Bangs of heavy ion collisions. 

2.1 Single particle spectra and two particle correlations 

The invariant single particle spectrum is obtained P, 0] as 


iVi(k) = 


cPn 
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27rmtdmt dy (2vr)'' 
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exp I --— I s 


(2.7) 
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The correlation function is found in the binary source formalism 

C 2 {ki,k 2 ) = 1 + A*H(Q||) exp - QIrI - QlR 


as: 


( 2 . 8 ) 


The pre-factor H((5||) of the BECF induces oscillations within the Gaussian envelope as 
a function of Qy. This oscillating pre-factor satisfies 0 < H(Q||) < 1 and H(0) = 1. In 
practice, the period of oscillations is larger than the corresponding Gaussian radius, so 
the oscillations are difficult to resolve. The above invariant BL-H form of the two-particle 
correlation function can be equivalently expressed in the frequently used but not invariant 
Bertsch-Pratt (BP) form in the LGMS frame [17|, within the H = 1 approximation; 


C’ 2 (ki, ka) = 1 + A* exp [-R^Q^ - RIqI - RjOf - 2RliQoQi] , 


(2.9) 


where the dependence of the fit parameters on the value of the mean momentum of the pair 
is suppressed. The above formulas for the BEGF and IMD, as were used in the fits, have 
been introduced in refs. and summarized recently in ref. |8|. We recommend 

this latter review paper for the formulas that relate the BL-H model parameters to the 
above forms for the spectra and correlation functions, in particular, eqs. (84-105), (115-118) 
and (129-140) of ref. [ 0 ]. 


3. Fitting preliminary STAR and PHENIX data on An -|- An at RHIC 

Here, we reconstruct the space-time picture of particle emission in Au -|- Au collisions 
at RHIG within the BL-H framework, by fitting simultaneously the PHENIX and STAR 
preliminary data on two-particle correlations and single-particle spectra as presented at 
the Quark Matter 2001 conference, refs. |jl^, H, | 


. For a proper core-halo correction 
oc 1/y/Xf the experimental values of the intercept parameter X:t,{y,mt) have to be taken 
from the measurements. In the lack of these values in ref. |^, we have utilized 

their average A* for a core-halo correction when fitting the PHENIX spectra. In particular, 


- 3 - 


PrHEP hep2001 














International Europhysics Conference on HEP 


Andras Ster 


the following average values were used for the various particle types: A*(7r) = 0.39 ± 0.14, 
Xif{K) = 0.80 (estimated from the NA44 data on kaon-kaon correlations at CERN SPS 
A* (p) = 1 (the fraction of long lived resonances that decay to anti-protons is neglected). In 
case of the STAR data, we have utilized the same values for A*(AT) and A*(p) however, for 
pions we have utilized the A*(mt) values of ref. [^. Note also that we have performed the 
data analysis within the = 1 approximation. Here we improve on our earlier results |25| 
by taking into account an mt dependent core-halo correction for the STAR spectra and 
correlations, and by fitting the absolute normalization of the single particle spectra in both 
experiments, properly utilizing the fugacity and quantum statistical factors. This allows 
us to extract the chemical potencial in the center of the fireball, in contrast to our earlier 
hts 1^1 where the absolute normalization of the particle spectra and the central value of 
the chemical potential distribution were not yet determined. Unique minima are found 
and a good is obtained for both data sets. Within errors, all the fit parameters 

remained the same as in ref. 1^ , but x^/NDF decreased slightly. See Figs. 1 and 2 for 
an illustation. The hypothesis that pions, kaons and protons are emitted from the same 
hydrodynamical source is in a good agreement with all the fitted data. 


4. Conclusions 

We find that the PHENIX and STAR data on single particle spectra of identified 7r“, K~ 
and p as well as detailed mt dependent HBT radius parameters are consistent with the 
Buda-Lund hydro model as well as with one another. The final state of central Au -|- 
Au collisions at RHIC corresponds to a cylindrically symmetric, large {Rg = 7.3 ± 0.7 
fm) and homogenous (Tq = 142 ± 4 MeV) fireball. A large mean freeze-out time, tq = 
8.4 ± 0.5 is found with a short duration of particle emission and small inhomogeneities of 





BL-Hydro 

STAR 

PHENIX 

Au-bAu 

Pb-bPb 

h+p 

parameters 

preliminary 

preliminary 

{RHIC) 

{SPS) 

SPS 

To [MeV] 

144 

± 5 

139 

± 5 

142 

± 4 

139 

± 6 

140 

± 3 

(ut) 

0.86 

± 0.10 

0.68 

± 0.3 

0.77 

± 0.09 

0.55 

± 0.06 

0.20 

± 0.07 

Rg [fm] 

8.0 

± 0.5 

6.6 

± 0.3 

7.3 

± 0.7 

7.1 

± 0.2 

0.88 

± 0.13 

To [fm/c] 

8.9 

± 0.5 

7.9 

± 0.3 

8.4 

± 0.4 

5.9 

± 0.6 

1.4 

± 0.1 

At [fm/c] 

0.5 

± 1.0 

0.6 

± 1.2 

0.5 

± 1.1 

1.6 

± 1.5 

> 1.3 

± 0.3 

Ay 

1.0 

± 0.1 

1.5 

± 0.1 

1.2 

± 0.3 

2.1 

± 0.4 

1.36 

± 0.02 

{^)r 

0.09 

± 0.01 

0.04 

± 0.01 

0.06 

± 0.03 

0.06 

± 0.05 

0.71 

± 0.14 

{^)t 

1.6 

± 0.4 

0.86 

± 0.09 

1.2 

± 0.4 

0.59 

± 0.38 

- 


yf [MeV] 

0 

(fixed) 

0 

(fixed) 

0 

(fixed) 

- 


- 


[MeV] 

46 

± 11 

- 


- 


- 


- 


hS [MeV] 

300 

± 18 

376 

± 38 

338 

± 28 

- 


- 


xVndf 

32/54 

= 0.59 

46/58 

= 0.79 

0.69 

1.20 

0.94 


Table 1: Preliminary source parameters from simultaneous fittings of preliminary RHIC data of 
PHENIX and STAR on particle spectra and HBT radius parameters with the Buda-Lund hydrody¬ 
namical model. The third column indicates their average. The (average) fit parameters are shown 
for Pb-I-Pb collisions at CERN SPS Q and for h-bp collisions at CERN SPS [ p^ . 
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Figure 1: Simultaneous fits to STAR parti- Figure 2: Simultaneous fits to PHENIX 
cle spectra and HBT radius parameters. particle spectra and HBT radius parameters. 


the temperature profile during the particle emission process. At the reconstructed final 
state, the hadronic fireball expands three-dimensionally with a strong transverse flow, 
characterized by the Hubble constant H = {ut) = 0.77 ± 0.09. The major difference 
between the hnal state of heavy ion collisions at RHIC and at CERN SPS seems to be an 
increased freeze-out time and an increased transverse flow or Hubble constant at RHIC. 
Within the presently large errors we do not find a significant increase in the reconstructed 
geometrical source size when increasing the energy from SPS to RHIC. This and other 
questions are being addressed by an attempt to include the pseudo-rapidity distribution 
at RHIC into the fitted data sets, and by repeating the analysis using final, published 
PHENIX and STAR data. 


We have found a non-vanishing chemical potential for kaons and anti-protons in the 
center of the fireball from the absolutely normalized single-particle spectra, while the pion 
data were well described in all cases with a vanishing pion chemical potential in the center 
of the fireball. These values together with the inhomogeneous chemical potential distribu¬ 
tion of eq. (2.4) indicate a clear deviation from chemical equilibrium in the reconstructed 
hadronic final state. 

The similarities and the differences between an effective Quark Matter (QM) stage 
and a Quark Gluon Plasma (QGP) phase have been summarized recently in ref. 

The observed short duration of particle emission and the large transverse flow at RHIC 
contradicts to the picture of a soft, long-lived, evaporative Quark Gluon Plasma phase, that 
would consist of massless quarks and gluons. However, the final state does not exclude a 
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transient, explosive, suddenly hadronizing Quark Matter phase, that could be characterized 
by massive valence quarks, the lack of gluons as effective degrees of freedom, and a hard 
equation of state. 

This research has been supported by the Hungarian OTKA T026435 and T034269, the 
NWO - OTKA grant N025186, by a NATO Science Fellowship, by an US NSF - Hungarian 
MTA-OTKA grant and by the DOE grant DE - EG02 - 93ER40764. 
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